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Advanced Monte Carlo simulations are used to study the effect of nano-slit confinement on metric and topo- 
logical properties of model DNA chains. We consider both linear and circularised chains with contour lengths 
in the 1.2-4.8 /xm range and slits widths spanning continuously the 50-1250nm range. The metric scaling pre- 
dicted by de Gennes' blob model is shown to hold for both linear and circularised DNA up to the strongest 
levels of confinement. More notably, the topological properties of the circularised DNA molecules have two 
major differences compared to three-dimensional confinement. First, the overall knotting probability is non- 
monotonic for increasing confinement and can be largely enhanced or suppressed compared to the bulk case by 
simply varying the slit width. Secondly, the knot population consists of knots that are far simpler than for three- 
dimensional confinement. The results suggest that nano-slits could be used in nano-fluidic setups to produce 
DNA rings having simple topologies (including the unknot) or to separate heterogeneous ensembles of DNA 
rings by knot type. 



I. INTRODUCTION 

Nano-devices are increasingly used to investigate the 
impact of spatial restraints on the statics and dynamics 
of polymers and biopolymers. DNA is ideally suited 
for such studies both for its applicative potential (nano- 
confinement can be used to sort, sieve and sequence 
DNA 0J[4]l) and because its elastic properties can be ex- 
ploited to elucidate issues of primary interest in polymer 
science I5UTT1I. In fact, the width of presently-available 
confining nano-devices (a channel, a slit, a pit etc.) can 
be set to compete with one or more of the characteristic 
lengths of a DNA molecule and hence allow for probing 
different physical regimes l8l fTTHT8l . 

A classic example is offered by the one-dimensional 
(channel) confinement of linear DNA molecules whose 
elongation follows two different scaling laws, named af- 
ter de Gennes [12] and Odijk (T9j, according to whether 
the channel width is larger or smaller than l p l20l . A 
second notable instance is offered by DNA molecules 
confined in three-dimensions, as inside a viral capsid. 
In the latter case, the interplay of the genome length 
and capsid size, can profoundly affect the entanglement 
of the packaged DNA l2Bi24ll . Both aspects are aptly 
revealed by the abundance and complexity of knots that 
are trapped in DNA filaments that circularise inside vi- 
ral capsids lfTTl[T7ll25H30ll. 

Compared to the above one- and three-dimensional 
cases, two-dimensional polymer confinement is much 
less characterised both for the metric and topological 
properties. To the best of our knowledge there exist only 
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a limited number of experimental and numerical stud- 
ied of the size and shape of DNA molecules in nano- 
slits [l8l l3Tti33ll . These seminal studies have ad- 
dressed the interesting issue of whether multiple scal- 
ing regimes exist in linear DNA filaments inside slits. 
It appears that the same question has not been posed 
for circularised DNA molecules where, unlike the linear 
case, entanglement is trapped in the form of topological 
constraints (knots) which affect the physical behaviour 
and biological functionality (34ti4T1l . The impact of slit 
confinement on the topology of circularised polymers 
has been previously addressed by early simulations of 
lattice polymer models [42 ] so that its implications for 
slit-confined DNA are virtually unexplored. 

As a step towards clarifying the above-mentioned as- 
pects we carry out a systematic numerical study of lin- 
ear and circularised chains in nano-slits whose salient 
physical properties, namely the thickness and bending 
rigidity are set to match those of dsDNA. Specifically 
we use advanced Monte Carlo techniques to character- 
ize the shape, orientation and entanglement of linear 
and circularised chains of contour lengths in the 1.20- 
4.8 fim range and confined in nano-slits of width from 
50 to 1250 nm. 

The computational scheme is first used to charac- 
terize the metric scaling by varying independently the 
chain contour length and the slit width. The metric 
scaling predicted by de Gennes' blob model is found 
to hold, up to the strongest levels of confinement, for 
both linear and circular chains. 

Next, we show that two-dimensional confinement 
affects the ring topology in a qualitatively different 
way from three-dimensional confinement. In fact, 
for increasing three-dimensional confinement both the 
incidence and complexity of knots grows monotoni- 
cally (23J [24). In the two-dimensional case, instead, 
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the knotting probability is non-monotonic and the knot 
spectrum mostly consist of the simplest knot types. The 
fraction of knotted molecules is maximum, and exceeds 
by several times the bulk (unconstrained) value, when 
the slit width is a fraction of the average extension of the 
unconstrained chains. A reduction of the width below 
this value, causes the knotting probability to decrease 
rapidly. 

This fact offers the opportunity to use the set-up of 
two-dimensional confinement in nano- slits as an effec- 
tive method to produce DNA rings with an a priori con- 
trol on their topology and that can be tuned to yield a 
proportion of knotted molecules that is much higher or 
much lower than the bulk one. In addition, the sensitive 
dependence of the conformational entropy of knotted 
molecules with different topology on the size of the con- 
fining region, suggests that width-modulated slits could 
be profitably used in nano-fluidic setups to sort circular 
DNAs by knot type. 



II. METHODS 

A. The model 

For this study, dsDNA is modelled as a semi-flexible 
chain of cylinders. This coarse-grained model cap- 
tures appropriately the metric and topological proper- 
ties of DNA both when it is unconstrained (43J 04) 
and when it is packed at moderate densities, that is 
for interstrand separations larger than « 5nm l45l . In 
such conditions, which are met in our study, dehydra- 
tion effects and cholesteric interactions can be neglected 
El El EH SSI S3- It is further assumed that DNA is 
in a concentrated buffer of monovalent counterions so 
that the screened electrostatic repulsion can also be ne- 
glected. The cylinder diameter is accordingly set equal 
to the dsDNA hydration diameter, d = 2.5nm. No- 
tice that no potential energy term is introduced to ac- 
count for the DNA torsional rigidity. Its interplay with 
bending rigidity can, a priori impact the relevant con- 
formational space of densely packed DNA, as argued in 
ref. 1 48 ] for closed DNA chains. However, recent in- 
dependent simulations of three-dimensionally confined 
DNA with free ends have shown that torsional effects 
are negligible even at packing densities higher than con- 
sidered here |49l . We shall therefore work under the hy- 
pothesis that the slit-confined DNA chains are torsion- 
ally relaxed both in the linear and in the circularised 
form. For the latter case it is assumed that DNA cir- 
cularization occurs via non-covalent annealing of the 
DNA sticky ends so that all torsional stress can still be 
released. 

We consider open and closed chains of contour length 
equal to L c ranging from 1.2 fim to 4.8 fim and describe 
them as chains of cylinders with N = 120 .. . 480 ver- 
tices, respectively, as shown in Table [I] This discretiza- 
tion level corresponds to a long axis of the cylinder, b, 



approximately equal to lOnm which is a fraction of the 
DNA persistence length, l p = 50nm. 



N 


bp 


L c [nm] 


Linear Rg [nm] 


Circular Rg [nm] 


120 


3534 


1201 


128 ±3 


94 ±3 


240 


7068 


2403 


190 ±5 


136 ±5 


320 


9426 


3204 


222 ±6 


159 ±6 


360 


10602 


3604 


236 ±6 


169 ±6 


400 


11782 


4006 


251 ±6 


179 ±6 


480 


14136 


4806 


276 ±7 


198 ±7 



TABLE I: Parameters of the simulated linear and circular 
DNA molecules: N is number of cylinders of the model 
chains, bp is the total number of base-pairs and L c is the cor- 
responding contour length. The last two columns provide the 
bulk values of the root mean square radius of gyration, Rg of 
the linear and circular chains. 

The potential energy of a chain configuration, T, con- 
sists of two terms, E r = E exc i_ vo i(Y) + E b (T). The 
first term is the excluded volume interaction which en- 
forces the chain self-avoidance. E exc \- vo \ is set equal 
to "infinity" if two non-consecutive cylinders overlap, 
otherwise is set equal to zero. The second term, E b is 
the bending potential energy: 

E b = -K b T l ^^2u-U +1 (1) 

i 

where ti is the orientation of the axis of the ith cylinder 
and the temperature T is set to 300K. For linear chains, 
i runs from 1 to N. For circular ones, instead, i runs 
from 1 to N + 1 with the proviso that t/v+i =t\. 

DNA confinement in nano-slit is modelled by setting 
an upper bound, D, to the calliper size of the chains 
along one direction, conventionally taken as the z axis 
of the system Cartesian coordinate, see Fig. [I] Because 
the z direction is perpendicular to the two impenetra- 
ble planes delimiting the slit, all observables measured 
along this direction will be indicated with the _L sub- 
script, while those measured parallel to the slit plane 
will be denoted with the 1 1 subscript. 

B. Advanced sampling techniques 

An advanced Monte Carlo sampling technique is 
used to compute the equilibrium metric and topolog- 
ical observables of open and closed chains subject to 
slit confinement. The technique is described in detail in 
a recent review article |11] and therefore we shall de- 
scribe it only concisely. 

The method consists of using a standard set of 
crankshaft and pivot moves to evolve the chains but 
the canonical weights used for the Metropolis accep- 
tance criterion are generalised to ensure that configu- 
rations are sampled according to the following (non- 
normalised) statistical weight: 

J E excl _ vol (T) + E b (T) - /iA±(T) \ ^ 
eXP { )• (2) 
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where Aj_(r), is the configuration calliper size mea- 
sured along the confining, z direction. The non-negative 
quantity /i can be viewed as an anisotropic pressure pro- 
moting the confinement along z. In fact, by increasing 
/i one can reduce the average calliper size of the gener- 
ated configurations. By using a multiple Markov chain 
scheme, consisting of several parallel Monte Carlo sim- 
ulations each at a different value of fi (23N50), it is pos- 
sible to sample very effectively configurations covering 
a wide range of calliper size, A^(T). Finally, a ther- 
modynamic reweighting technique analogous to the one 
used in ref. l23l is used to combine the observables com- 
puted for configurations generated at various values of 
ji. This procedure gives an estimate of the equilibrium 
properties of chains that can be accommodated inside 
slits of width, D, in the 50nm - 1250nm range. 

C. Metric properties 

The average size, shape anisotropy, and orientation 
of a given ring are characterized by means of the eigen- 
values and eigenvectors of its gyration tensor, Q. The 
entries of this 3x3 symmetric matrix are given by: 

1 N 

Qa,f3 = ^ ^2( r i,a ~ ra)(r i)/3 - Tp), (3) 

where is the position of the ith vertex of the chain, 
a and f3 run over the three Cartesian components and 
r = i*i is the average vertex position in the chain. 

The non-negative eigenvalues of Q, ranked for decreas- 
ing magnitude, are indicated as Ai, A 2 and A 3 and their 
associated (normalised) eigenvectors are indicated as 
e\, e*2 and e^. 

The sum of the eigenvalues provides the square ra- 
dius of gyration of the chain: 

R 2 g = Ai + A 2 + A 3 , (4) 



while the adimensional ratios A1/A3 and A2/A3 are 
used to characterize the chain shape anisotropy. 

The chain orientation inside the slit is conveniently 
described by the angle a formed by each principal gy- 
ration axis and the slit plane: 

cos ai = \/l - \el • z| 2 (5) 

where z is the unit vector in the z direction. 

To parallel the quantities that are customarily mea- 
sured in nano-slit experiments, for each chain configu- 
ration we compute its mean square size projected in and 
out of the slit plane which are, respectively: 

R\\ = ^2 ^i C °S 2a ii ( 6 ) 
i=l,2,3 

Rl = R 2 g - Rf\ = A *(! - cos2 a *) ' W 

z=l,2,3 

D. Scaling properties and the deGennes regime 

The de Gennes' blob theory provides a powerful 
reference framework for characterizing spatially con- 
strained chains lfT2l . The theory, originally formulated 
to improve on standard mean-field theories for semi- 
dilute polymer solutions and HUES), ^ as P rove d use- 
ful to understand the behaviour of polymers confined in 
regions with transverse size larger than the molecule's 
persistence length. 

When such condition is met, the theory of de Gennes 
proposes to describe the self-avoiding chain, formed by 
N monomers of size b, as a sequence of "blobs", with 
each blob consisting of g monomers and having linear 
size about equal to the slit width, D, see Fig. [2] Within 
each blob, the chain does not experience the confining 
constraints and hence it behaves as an unconstrained 
self-avoiding walk: D ~ g U3JD , where v^d ~ 0.588 
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FIG. 2: Pictorial representation of a confined polymer as a 
string of DeGennes' blobs 



is the three-dimensional self-avoiding walk scaling ex- 
ponent (53) . At the same time, the string of blobs can be 
viewed as being a two-dimensional self-avoiding walk, 
with the blobs being the effective monomeric units of 
the walk. Accordingly, the in-plane size of the chain ex- 
pectedly scales as R\\ ~ (N/g) V2 ^D where v?d =3/4 
is the self-avoiding exponent in two dimensions. 

This relationship can be recasted as R\\ ~ 
(N/g) v ™D ~ A/^2d£)-^dMd£) or? equivalent^, 

R\\ /^X 1 — W-3D , d k -0.276 

where we have exploited the proportionality of R° and 
iijj which are respectively the root mean square values 
of the radius of gyration and the in-plane projected size 
of unconstrained chains. It is worth pointing out that 
the pref actor entering the scaling relationship of eq. [8] 
changes according to whether the blobs are isometric 
or not (13] [20). This condition depends, in turn, on 
whether the slit width D is larger or smaller than / d. 
For our case, the latter quantity is equal to Ifim and 
therefore is always larger than the widths of the slits 
considered. The blobs are therefore expected to be ani- 
sometric (extended de Gennes regime) (20). 



E. Topological properties 

The characterization of the entanglement of circular 
configurations is completed by establishing their knot- 
ted state. It should be noted, in fact, that the crankshaft 
moves employed for the MC evolution preserve the ring 
connectivity but not necessarily its topology [54]. In 
fact, even if the initial and final state of a crankshaft 
move are self-avoiding, the virtual crankshaft move- 
ment that bridges them can entail several bond crossings 
and hence can change the ring knotted state. 

The topology of a ring is established using the same 
methodology of ref . (23] |45) which is articulated over 
the following steps. First, the ring geometry is sim- 
plified by topology preserving moves fTT1l23ll55ll56) . 
These moves consist of ring deformations involving no 
bond crossing and promote the collinearity of triplets 
of ring vertices that are nearby in sequence. After the 
chain simplification, a two-dimensional projection of 



the resulting configuration is next taken, encoded as 
a digital Dowker code, simplified algebraically and fi- 
nally compared against a lookup table of prime knots 
with up to 16 crossings. If a match is found for all the 
prime component of a knot, then the ring knotted topol- 
ogy is unambigously established otherwise it is labelled 
as "unknown". 

Notice that the "unknown" states may correspond to 
complicated knots with prime components having more 
than 16 crossings, or may be a simple knot that could 
not be simplified enough to positively identify it. 

Further insights on the interplay between the topo- 
logical properties of circular DNAs and the degree of 
confinement may be obtained by considering the av- 
erage length of the knotted region in the ring. Locat- 
ing the knotted portions in closed chains is a challeng- 
ing problem, that presents the same conceptual difficul- 
ties as establishing the knotted state of a geometrically- 
entangled open chain Bj][54). Although the applica- 
tion of topological concepts to open chains is a priori 
ambiguous, it has been recently verified that very differ- 
ent methods for locating knots in linear chains yield re- 
markably consistent results [57 ]. Among them we men- 
tion those based on the statistical closure scheme (58) 
and the minimally-interfering one (57) which is used 
here for its robustness and numerical efficiency. 



III. LINEAR DNA IN NANO-SLITS 
A. Effect of confinement on chain orientation and shape 

The effect of slit confinement on the chain orientation 
is illustrated in Fig. [3^. 

The data reflect the in-plane and out-of-plane orien- 
tations, see eq. [5j of the three axes of inertia of open 
chains of contour length, L c = 3604 nm. The bulk 
root mean square radius of gyration of these chains is 
Rg = 236 nm (see Table [i]) and its location in the plots 
of Fig.[3]is depicted by the vertical dot-dashed line. For 
slit widths much larger than the chain bulk size the three 
mean square cosines are all close to the value 1/3, as 
expected in case of isotropic orientation of the chain in 
space. 

The orientation of the longest axis of inertia begins 
to be noticeably non-isotropic when the slit width is 
about 700nm, which exceeds three times the bulk ra- 
dius of gyration. This is consistent with the strong shape 
anisotropy of unconstrained chains (A1/A3 = 13.8 and 
A2/A3 = 3.05) which implies that the longest axis of 
gyration cannot be freely-accommodated inside the slit 
already when the width D is several times larger than 
R? g . The second longest axis of gyration is oriented 
isotropically, and maintains its bulk value, down to slit 
widths D « 2R? g = 472nm (shown in the panels of 
Fig. [3] as a dashed vertical line). The value of A2 re- 
duces appreciably when D is decreased down to about 
Rg and remains approximately constant for smaller slit 
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widths. Analogous properties are found for the other 
chain lengths considered in this study. 

For D < Rg, both the first and the second axes of 
inertia are mostly confined in the slit plane, while the 
third is oriented perpendicularly to it. 

The progressive in-plane orientation of the two prin- 
cipal axes of inertia is paralleled by the flattening of the 
chain in the slit plane, see Fig. [3]). From Fig. |3j: it is 
noticed that for D < 2R g the progressive in-plane ori- 
entation of the second axis of inertia is accompanied by 
a decrease of its eigenvalue A 2 to a limiting value that 
is smaller than in the unconstrained, bulk case. The op- 
posite is true for Ai. As a result, strongly confined open 
chains tend to have an elongated shape, as in the exam- 
ple of Fig.[lJ. 

The adimensional data in Fig. [3}),c compares 
favourably with available numerical and experimen- 
tal results for the limiting cases of no- and strong 
two-dimensional confinement. In particular, Haber et 
al. l59l report the aspect ratios ^/Ai/As = 4.1 and 
Y^Ai/A2 = 2.3 for fully-flexible unconstrained poly- 
mers, which are close to the corresponding values 3.7 
and 2.1 found here. Finally, Maier et ah, in their study 
of A-phage DNA adsorbed on a flat bilayer lipid mem- 
brane E3, found A 1 /R 2 g = 0.84 and A 2 /R 2 g = 0.14 
which is in accord with the values A\/R 2 = 0.84 ± 
0.02 and A 2 /R 2 g = 0.148 ± 0.008 found here for the 
strongest two-dimensional confinement. 



B. Effect of confinement on chain size 





200 400 600 800 1000 1200 
D (nm) 

FIG. 3: Metric properties of a linear chain L c = 3604 nm 
inside a slit of width D. The orientation, shape and size of the 
gyration ellipsoid of the confined chain are shown in panels 
(a-c). The root mean square values of the chain radius of 
gyration, and its in- and out-of-plane projections are shown 
in panel (d). The range and scale of the x axis (representing 
the slit width, D) is common to all panels. The two vertical 
lines mark the reference slit widths R° g (dot-dashed) and 2R Q g 
(dashed). 



The global, in-plane and transverse size of the chain 
at various levels of confinement are shown in Fig.[3]i 

The transverse size, R± and the parallel one, R\\, 
have opposite monotonic dependence on D which re- 
sult in a non-monotonic behaviour of R g as a function 
of D. From Fig.[3]l, it is seen that R g is minimum for D 
about equal to the bulk radius of gyration of the chain, 
Rg. This property is found consistently across chains of 
all considered lengths, as shown in Fig. [4^. 

To test the validity of the de Gennes' blob regime we 
analyzed the D dependence of R\ | . The result is shown 
in Fig. |4j) which reveals a remarkable collapse of data 
points from curves of lengths L c between 2.4/im and 
4.8/im. As shown in the figure, for high confinement 
R\\ has a power law dependence on D, with exponent 
equal to —0.27 ± 0.05. This exponent, estimated by fit- 
ting the most asymptotic data i.e. L c > 3604 nm, is 
well compatible with the -0.276 value predicted by the 
de Gennes' theory. Notice that, at all considered lengths 
no deviation from the de Gennes' theory is observed, 
consistently with the recent experiments of Doyle et 
al. fT8) on longer DNA filaments. 
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can be directly compared with those in Fig. [3] for an 
equally-long open chain. 
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FIG. 4: (a) The minimum three-dimensional size of the chain 
is achieved when the slit width is about equal to the average 
radius of gyration of the unconstrained chain, (b) Below such 
slit width the relationship between the average in-plane size, 
R\\, chain contour length, N, and slit width, D, is well ac- 
counted for by the DeGennes blob theory 



IV. CIRCULAR DNA IN NANO-SLITS 

We consider DNA rings with the same contour length 
as the previously-discussed open chains, see Table [TJ 
and for the same range of confining slit widths. The 
properties of rings are first characterised at the met- 
ric level and then at the topological one. We shall 
particularly focus on the type, abundance and size of 
the various knots that are trapped in confined model 
DNA molecules that are circularized in equilibrium. 
These aspects are known to affect significantly the static 
and dynamics of polymers both in dilute (62 |62) and 
dense solutions 1 63 ] or tightly packed in small vol- 
umes O El E3 E3 ES S3 or pulled through narrow 
channels (27l |28) and are hence expected to be impor- 
tant for two-dimensional confinement too. 



A. Effect of confinement on ring orientation and size 

The dependence of a ring size and orientation on the 
confining slit width is illustrated in Fig. [5] The data 
pertains to rings of length L c = 3604 nm and therefore 
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FIG. 5: Metric properties of a circular chain L c = 3604 nm 
inside a slit of width D. The orientation, shape and size of the 
gyration ellipsoid of the confined ring are shown in panels (a- 
c). The root mean square values of the ring radius of gyration, 
and its in- and out-of-plane projections are shown in panel (d). 
The range and scale of the x axis (representing the slit width, 
D) is common to all panels. The two vertical lines mark the 
reference slit widths R° g (dot-dashed) and 2R° g (dashed). 
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The comparison of panels (b) and (c) in the two fig- 
ures indicates that the shape of closed chains is signif- 
icantly more spherical than open ones at all levels of 
confinement. The effect is aptly illustrated by the typi- 
cal open and closed chains shown in Fig.[T] 

Notwithstanding the different shape, it is found that 
the changes in the ring orientation and size occur for 
values of the normalized slit width, D / 'R® that are sim- 
ilar to the linear case. In particular, the data in panel 
(a) show that the in-plane orientation of the ring prin- 
cipal axis of gyration sets in for D « 2R g while for 
the second one occurs at D w R g . Furthermore the 
average ring size, R g , shown in panel (d) has a non- 
monotonic dependence on D and its minimum value 
(again given by the opposite trends of R\\ and R±) oc- 
curs for D w R° as for the linear case. 



B. Scaling properties and the deGennes regime 

Because the bulk metric scaling exponents v?d and 
^3D are the same for open and closed chains, it is ex- 
pected that the de Gennes scaling regime, originally 
based on the blob argument for linear chains, ought to 
hold for circular ones too. 

However, the closed character of the rings introduces 
correlations in the orientation in pairs of bonds at a large 
chemical (contour) distance that have no analog in lin- 
ear chains (64j [53 

To clarify the effect we show in Fig. [6] the same ob- 
servables reported in Fig. [4] to analyze the de Gennes 
regime for confined linear chains. 

Panel (a) shows that, except for the data at the small- 
est ring contour length, there is a good overlap of the 
rescaled ring size versus width curves. In particular, for 
all considered chain lengths the minimum overall ring 
size is obtained for D « R®, consistently with the open 
chain case 

For D < R° g the mean in-plane size shows a be- 
haviour that is compatible with the de Gennes scaling 
regime, as illustrated in panel (b). The power-law fit of 
the asymptotic data in the figure, i.e. L c > 3604nm and 
D/R° g < 1, yields the scaling exponent -0.26 ± 0.04. 
This value is compatible with the -0.276 value predicted 
by the de Gennes theory, consistently with what found 
for the linear case. 



C. Topological entanglement 

1. Knotting probability 

The Monte Carlo strategy employed here and in other 
previous studies provides the means to characterize the 
occurrence of non- trivial knots in rings that are circu- 
larised in equilibrium. The knotting probability of un- 
constrained rings depends on several lengthscales such 
as the molecule's persistence length, its contour length 



DC 



DC 



1.6 
1.4 

D 

1.2 



Circular DNA 



1- 



( a ) 0.8 



1.6 




1201 nm 
2403 nm 
3204 nm 
3604 nm 
4006 nm 
4806 nm 



D / R 



o _ 

DC" 
DC" 



(b) 



1.2 





— 1201 nm 




— 2403 nm 




3204 nm 




v — 3604 nm 




S^v — 4006 nm 




4806 nm 




♦^S. -0.26 






1 1 1 


1 1 1 1 1 



1 

D/R 



10 



9 



FIG. 6: (a) The minimum three-dimensional size of the ring is 
achieved when the slit width is about equal to the average ra- 
dius of gyration of the unconstrained ring, (b) Below such slit 
width the relationship between the average in-plane size, R\\, 
ring contour length, N, and slit width, D, is well accounted 
for by the de Gennes' blob theory. 



and thickness (66] I67I . Accordingly, the topological 
complexity of the circularised molecules can be tuned, 
to some extent, by varying these lengths, for example 
by intervening of the solution ionic strength |43ll68l . 

Spatial confinement introduces a further lengthscale 
in the problem, the width of the confining region, which 
can affect dramatically the molecule's entanglement. 
The effect has been addressed by several studies of 
isotropic three-dimensional confinement (XT] |23] [25j 
l45ll which established that progressive reduction of the 
volume resulted in an ever-increasing fraction of knot- 
ted molecules and of their geometrical and topological 
complexity l24l . 

By comparison, the case of two-dimensional confine- 
ment is virtually unexplored, except for lattice polymer 
models | 42, 69, 70 ]. Making progress in this problem is 
not only important to fill a conceptual gap but is espe- 
cially so in view of the dramatic qualitative differences 
that are expected a priori compared to the isotropic, 
three-dimensional, confinement. 

In fact, it should be considered that for extreme 
two dimensional confinement, i.e. for slit width 
only slightly larger than the chain intrinsic thickness, 
over/under-crossings of the chain are impossible and 
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knots cannot be formed. Therefore, unlike the isotropic 
case, progressive slit confinement cannot result in an 
ever increasing complexity of the ring topology and ge- 
ometry. 

To clarify this aspect we investigated the dependence 
of the knotting probability on D. The results are shown 
in Fig. [7^ and indicate that the ring knotting probabil- 
ity is non-monotonic and for a certain slit width, D e 
reaches a maximum enhancement compared to the bulk 
case. Notice that the peak value exceeds the bulk knot- 
ting probability by several factors. The figure indicates 
that the limiting value of zero knotting probability is 
rapidly reached for D < D e . 
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FIG. 7: (a) Knotting probability of equilibrated rings of dif- 
ferent contour length, (b) Dependence on the slit width, D, 
of the average radius of gyration and knotting probability for 
a circular DNA ring of L c — 4806nm. The vertical dashed 
line marks the slit width, D* , corresponding to the minimum 
three-dimensional size of the ring. The dot-dashed line marks 
the maximum of the knotting probability. 



It could be intuitively expected that the maximum 
knotting enhancement occurs for the same slit width, 
D*, associated to the maximum chain density (mini- 
mum R g ) [42]. This is, however, not the case. In fact, in 
the example of Fig. [T]} the peak of the probability distri- 
bution occurs when the slit width is D e « 90nm, which 
is about one third of D*. As a matter of fact, over the 
considered range of contour lengths, D e is consistently 
equal to about half the bulk ring size. 



2. Knot spectrum 

It is most interesting to examine how the various, 
topologically inequivalent knot types, contribute to the 
nonmonotonic trend of the overall knotting probability. 

The results are shown in Fig. [8^,b which portrays 
the probability of occurrence of various knot types (in- 
cluding the unknot) for four different slit widths and 
L c = 4.8/im. It is seen that, at all level of confinement, 
the knot population consists of very simple knots, and is 
dominated by the simplest of all, the trefoil or 3i knot, 
whose peak abundance is ~14%. 

The result is notable in several respects. First, the 
knot spectrum of Fig. [8^ is qualitatively very different 
from the one of an equivalent model of circular DNA 
subject to isotropic, spherical confinement l45l . 
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FIG. 8: (a,b) Slit- width dependence of the knot spectrum of 
confined circular DNA with L c — 4806nm. D* is the slit 
width at which the mean 3D extension of the rings reaches its 
minimum value, (c) Maximun probability of occurrence for 
knot types 3i, 4i and 5i as a function of the contour length 
L c (dashed curves). Solid curves refer to the corresponding 
case of unconstrained circular DNAs. 

In fact, even for moderate three-dimensional confine- 
ment in a sphere of diameter equal to 1 80nm, the knot 
spectrum of molecules of Lc = 3.4/im is largely con- 
stituted by highly complex knots: trefoils are less than 
10 % of the non-trivial knots (which, in turn are 55% of 
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the ring population). Conversely, for slit confinement, 
at all considered values of L c and D, complex knots 
are virtually absent and the population of simple knot 
types, such as the trefoil, reaches values that cannot be 
matched by rings of the same or longer contour lengths 
subject to any degree of isotropic confinement (23). 

We finally turn to the dependence of the peak of the 
probability of occurrence of simple knots as a function 
of chain length. The data in Fig. [SJ; indicate a mono- 
tonic increase of the peak probability with L c . There- 
fore, the population of simple knots can be further in- 
creased compared to the values in panel (b) by suitably 
confining rings with L c > 5/im. At the same time, 
it is noticed that the derivative of the peak probability 
curves decreases with L c so that, by analogy with the 
bulk probability case, one can envisage that the curves 
will eventually decrease for sufficiently large values of 
L c . This is an interesting point that would deserve to be 
addressed in future studies. 



3. Knot length 

We conclude the topological characterization by re- 
porting on the degree of localization (in space and along 
the chain) of the knots formed in the confined rings. The 
knot location along the ring was established using the 
bottom-up search of the smallest arc(s) that, after clo- 
sure, has the same topology of the whole ring l57ll7TTL 
see the example in Fig. [9] 

The average contour length spanned by the knots (all 
non-trivial topologies) as a function of D is shown in 
Fig.flOlfor rings of length Lc = 4.8/mi. It is seen that 



up to the point where the average ring size is minimum, 
the length of the knotted region is about constant, and 
drops significantly at smaller slit widths. 

The same figure also shows the average in-plane 
and overall size of the knotted region, which is indi- 
cated with Rg and , respectively. It is apparent that 
the drop in knot length occurs when the slit width is 
about equal to the bulk value of the overall knot size, 
D & 125nm. These results suggest that for knotted 
polymer rings confined into slits the typical size of the 
knotted region introduces an additional length scale into 
the problem that could play a relevant role in under- 
standing the metric and mechanical properties of these 
confined molecules at equilibrium l70ll72l . 



V. CONCLUSIONS 

The equilibrium properties of linear and circularised 
model DNA molecules confined in nano- slits were stud- 
ied by numerical methods. By using thermodynamic 
reweighting techniques the properties of these sys- 
tems were systematically profiled as a function of the 
molecules' contour length and nano- slit width. The 
analysis proceeded by first considering the molecules' 



geometrical properties and next the topological ones. 
The former included the size, shape and spatial orien- 
tation of the molecules in the slit, while the latter con- 
sisted of the knotting probability and knot spectrum. 

We found that the scaling of the metric properties was 
described very well by de Gennes' blob theory for both 
linear and circular molecules of all considered lengths 
and down to the highest confinement. To the best of 
our knowledge, it is the first time that the applicability 
of the blob theory is discussed (and shown to hold) for 
equilibrated ring polymers. In addition, even for linear 
molecules, the fact that the de Gennes' scaling is inves- 
tigated upon varying both the molecules' contour length 
and the slit width represents a significant broadening of 
the scope of previous investigations. Our results are 
consistent with the conclusions of recent experiments 
of nano- slit confinement of linear A -DNA which found 
no indications of deviations from de Gennes' blob the- 
ory d. 

Finally, and most importantly, we characterized the 
topological entanglement that is trapped in molecules 
that circularise in equilibrium inside nano-slits. For 
both applicative purposes and for fundamental poly- 
mer science, the characterization of the physical prop- 
erties of DNA molecules with various topologies (lin- 
ear, circular planar, circular knotted) has attracted much 
interest in the past decades. Our contribution repre- 
sents a step forward in the characterization of such 
systems in the presence of confinement: while the 
three-dimensional confinement has been heavily inves- 
tigated (23] [25] [26] EU S3, the two-dimensional one 
by means of slits is still largely unexplored. Our study 
has pointed out a remarkable feature that has no parallel 
with three-dimensional confinement, namely the non- 
monotonicity of the knotting probability for increasing 
confinement. Specifically, it is seen that, among the 
rings resulting from equilibrium circularization (studied 
without enforcing topology conservation) the fraction 
of knotted rings has a maximum at a length-dependent 
slit width, D e . The peak knotting probability exceeds 
by several factors the one in the bulk. Upon decreasing 
D below D e , the knotting probability decreases very 
rapidly to zero. 

The above result shows that nano- slit confinement 
ought to be profitably used in contexts where an a pri- 
ori control of the knotted topology of DNA rings is de- 
sired. In particular, the non-monotonic dependence of 
the knotting probability as a function of confinement, 
makes it possible to conceive experiments where the slit 
width is tuned so that the knotting probability is either 
significantly enhanced or largely depressed with respect 
to the bulk case. 

The natural extensions of the present study are 
twofold. On the one hand it would be interesting to 
consider much longer chain lengths and establish the 
detailed length-dependence of D e and the associated 
maximum knotting probablity. On the other hand, it 
would be most interesting to repeat the present inves- 
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tigation to the one-dimensional confinement case (i.e. 
nano-channels) and analyze how the expected consecu- 
tive onsets of the de Gennes and Odijk regimes impact 
both the geometrical and the topological properties of 
confined rings. 



Reprinted with permission from 
C. Micheletti and E. Orlandini, Macromolecules, 2012, 
vol. 45, pp. 8668-8680. 
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FIG. 9: (a) Trefoil-knotted ring confined in a nanoslit. The configuration is the same as the one shown in Fig. |TJ;,d. The ring 
contour length is 0.48 /im and the slit width is 70nm. The shortest knotted portion of the ring is highlighted in red. (b) Close-up 
of the knotted portion; for graphical clarity the slit top is not shown. 



225 



E 




,i , 1 . 1 1 ' i i i i 1 

200 400 600 800 1000 



D (nm) 

FIG. 10: Three-dimensional and parallel mean estension 
of the knotted portion of circular trefoil knots with L c = 
4806nm. Circles refer to the average knot size expressed in 
unit of cylinders. 



